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Summary- 

Response  of  a  metallic  capsule  to  a  transient  internal  pressure  is 
analyzed  adopting  two  models.  The  first  is  a  two-dimensional  (2-D)  model 
with  finite  extent  in  the  plane  of  the  capsule's  cross-section  and  infinite 
length  normal  to  that  plane.  The  second  is  a  three-dimensional  (3-D) 
axisymmetric  model  consisting  of  a  short  thin  cylinder  capped  by  thin  disks. 
In  both  models,  transfer  matrices  relate  state  vectors  of  tractions  and 
displacements  at  two  segment  boundaries.  The  time  interval  from  pressure 
initiation  to  failure  of  its  material  measures  the  capsule's  effectiveness 
as  a  confining  boundary. 
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1.  Introduction 

Ceramic  materials  are  used  to  harden  vehicles  and  structures  against 
ballistic  events.  These  materials  benefit  from  a  higher  compressive  strength 
compared  to  metals  yet  their  tensile  strength  is  reduced  by  brittle  failure. 
Following  impact,  a  ceramic  material  shocks  the  projectile  initiating 
damage.  As  time  progresses  and  tensile  release  waves  form,  micro-cracks 
propagate  and  multiply  in  the  material  culminating  in  total  failure.  Strain 
energy  in  the  virgin  material  then  converts  into  kinetic  energy  ejecting  the 
debris  away  from  the  projectile's  path.  Prolonging  the  time  interval  when 
comminuted  material  remains  in  this  path  after  failure  is  one  way  to  extend 
the  material's  utility  in  eroding  the  projectile. 

In  what  follows  a  quantitative  account  is  given  of  an  experimental 
project  duplicated  in  two  laboratories  and  aimed  at  evaluating  ballistic 
performance  of  stacked  configurations.  This  account  is  interesting  as 
unplanned  differences  in  experimental  setups  revealed  the  strong  effect  of 
confinement  on  penetration. 

Controlled  ballistic  experiments  on  periodic  stacks  of  ceramic  layers 
bonded  by  thin  polymer  films  were  performed  at  two  laboratories  named  "A" 
and  "B"  for  shortness  (Weber  and  Hohler  (1997),  Weber,  et .  al .  (1999)). 

Ballistic  performance  was  measured  by  the  residual  penetration  pr  of  the 

projectile  into  a  thick  Aluminum  backing  block.  In  this  way  pr  is  a  measure 
of  the  projectile's  residual  kinetic  energy.  All  experiments  consisted  of 
stacks  made  of  100mm  x  100mm  square  Aluminum  Nitride  (AlN)  tiles  and  38mm 
total  thickness  of  either:  1  x  38mm  monolith,  2  x  19mm  layers,  3  x  12.7mm 
layers,  or  6  x  6.35mm  layers.  All  layers  were  bonded  by  0.254mm  uniform 
polyurethane  film,  adhering  to  the  tiles  by  heating. 

In  Series  3c,  the  stacks  in  the  "B"  setup  were  confined  laterally  by  14cm  x 
14cm  square  steel  tubing  12.5mm  thick,  and  axially  by  a  5mm  thick  steel 
cover  plate  perforated  centrally  with  a  30mm  diameter  hole  allowing  the 
projectile  to  strike  the  stack  (Fig.  1(a)).  The  stack  rested  on  a  60mm 

Aluminum  backing  block  to  serve  as  witness  plate  for  pr  measurements.  The 

"A"  setup  was  laterally  confined  was  achieved  by  welding  four  6.35mm  plates 
along  the  edges  to  form  a  box  keeping  a  5mm  clearance  around  the  stack,  and 
covered  by  a  6.35mm  plate.  The  cylindrical  projectile  was  made  of  Tungsten 
alloy  with  8.33mm  diameter  and  length  to  diameter  ratio  of  6,  weighing  50g, 

with  velocity  near  1170  m/s.  Fig.  1(b)  plots  pr  versus  number  of  layers  in 
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the  stack  for  all  "A"  and  "B"  stacks  tested.  In  this  Figure,  data  from  each 
series  of  experiments  is  labled  by  a  different  symbol  like  •,  A,  □,  O  and 
X.  At  first,  results  from  "A"  and  "B"  did  not  match.  "A"  measured  pr= 30mm 
for  the  1.5"  monolith,  and  pr  =0mm  for  the  1/2"  layer  stack  (  •  Fig.  1(b)). 

«B"  measured  /?r=0mm  for  the  38mm  monolith  and  6mm  for  the  12.7mm  layer 
stack  (not  shown  in  Fig.  1(b)).  A  review  of  the  "A"  experimental  setup 
revealed  that  lateral  confinement  failed  at  the  corner  welds,  ejecting  the 
side  plates  laterally  immediately  following  impact  thus  negating  the  effect 
of  confinement.  This  confirmed  that  the  principal  cause  for  the 

substantially  lower  pr  in  the  "B"  setup  was  the  strong  confinement  that 

delayed  removal  of-  comminuted  material,  from  the  projectile's  path, . 

increasing  its  erosion.  This  effect  masked  the  real  ballistic  performance  of 
ceramic . 

The  "A"  and  "B"  setups  repeated  Series  3c  experiments  with  the 
following  changes.  The  size  of  the  confining  tube  in  "B"  was  increased 
allowing  a  5  mm  clearance  around  the  stack  for  lateral  expansion  of  the 
damaged  ceramic  without  interference  from  the  confining  tube.  Also,  a  3mm 
gap  was  left  between  cover  plate  and  top  face  of  the  stack.  In  the  "A" 
setup,  the  gap  between  stack,  sidewalls  and  cover  plate  was  increased  to 
20mm  practically  eliminating  the  effect  of  confinement.  Indeed,  results  for 

the  12.7mm  layer  stack  agreed:  "B"  measured  pr=  20mm  and  23mm  (A  Fig.  1(b))  , 
and  "A"  measured  /3r=20mm  and  26mm  (  •  Fig.  1(b)).  For  the  38mm  monolith  the 

mismatch  in  pr  persisted:  "B"  measured  pr  =6 . 8mm  and  Omm  (  A  )  pointing  to  a 

possible  problem  in  the  "B"  setup.  "B"  was  repeated  in  3  more  tests  on  the 

38mm  monolith  yielding  pr=  11  mm,  33mm  and  46mm  (  O  Fig.  1(b)).  This  unusual 

scatter  among  "B"  results  for  the  38mm  monolith  confirmed  that  a  fundamental 
difference  existed  between  the  two  setups  meriting  a  closer  look  at  the  "B" 
setup.  Examining  the  "B"  cover  plates  of  the  last  3  tests  revealed  that  the 

cover  plate  corresponding  to  pr  =llmm  suffered  substantial  plastic 
deformation,  while  that  corresponding  to  pr= 46mm  had  very  little  plastic 
deformation.  This  suggested  that  the  large  scatter  in  pr  was  caused  by  how 
the  cover  plate  confines  ejecta  in  the  direction  of  impact.  The  3mm  gap 
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between  cover  plate  and  stack  in  "B"  was  small  enough  to  confine  ejecta  at 
times  when  they  did  not  escape  through  the  central  hole  in  the  cover  plate 
remaining  in  the  path  of  the  projectile.  This  increased  resistance,  reducing 

p  and  creating  sufficient  pressure  to  cause  the  cover  plate  to  yield.  When 
ejecta  escaped  through  the  hole,  the  projectile  encountered  less  resistance 
increasing  p  and  keeping  the  cover  plate  intact.  This  demonstrated  that 

even  the  weak  confinement  of  the  cover  plate  randomly  affected  ballistic 
performance  by  retaining  or  allowing  ejecta  to  pass  through  the  hole. 

An  explanation  was  still  needed  to  account  for  the  match  in  results 
and  reduced  scatter  in  the  "B"  data  of  Series  3c  for  the  multiple  layer 
configurations  in  spite  of  the  cover  plate's  closeness.  That  explanation 

lies  in  the  concept  of  phase  velocity  C  .  For  the  monolith,  the  rate  of 

formation  of  comminuted  material  in  ejecta  is  approximately  the  compression 
speed  of  sound,  which  for  AlN  is  lOkm/s.  For  a  stack  made  of  12.7mm  ceramic 

layers  bonded  by  0.254mm  polyurethane  film,  C p  ~-<J(Ebhc )  /  ( pj\ )  =  1 . 7km/s 

where  Pc,hc  are  ceramic  mass  density  and  layer  thickness,  and  Eb,hb  are  film 

modulus  and  thickness  (El-Raheb  (1997)).  The  reduction  in  C  caused  by 

dispersion  from  weak  coupling  reduces  the  rate  of  ejecta  formation  by  a 
factor  of  6.  This  raises  the  probability  of  ejecta  escaping  from  the  hole 

and  in  turn  diminishes  scatter.  In  fact,  the  reduction  in  pr  scatter  for 
stacks  with  thinner  ceramic  layers  is  consistent  with  the  corresponding 
reduction  in  C  . 

To  bring  the  two  setups  into  equivalence  required  adjusting  the  "B" 
setup  in  Series  3d  by  increasing  the  distance  between  cover  plate  and  stack 
from  3mm  to  50mm.  This  allowed  sufficient  space  for  the  ejecta  to  expand 
without  restriction  from  confinement.  With  this  final  setup,  experiments 
measured  a  higher  pr  in  two  tests  with  the  38mm  monolith:  29mm  and  41mm 

(  □  Fig.  1(b)) .  Faced  with  this  unexpectedly  strong  effect  of  confinement, 

the  confining  tubing  in  the  "B"  setup  was  increased  from  14cm  to  19.5cm  and 
2  more  tests  were  performed  for  each  stack  reducing  the  scatter  to  less  than 

15%  in  Series  3e.  For  the  38mm  monolith,  pr=  25mm  and  31mm  (X  Fig.  1(b))  . 
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Note  the  large  scatter  in  pr  for  the  monolith  when  gap  between  cover 
plate  and  stack  was  3mm  (  A  and  0  Fig.  1(b)),  and  the  reduced  scatter  of 

all  data  when  that  gap  was  increased  to  50mm  (  D  and  X  Fig.  1(b)) .  The 

scatter  diminishes  as  layer  thickness  in  the  stack  is  reduced. 

Examining  the  experimental  results  presented  above,  it  is  evident  that 
metallic  encapsulation  of  ceramic  tiles  may  extend  the  time  interval  of 
active  resistance  against  the  projectile  also  known  as  dwell  time.  One 
drawback  is  added  weight  implying  that  an  optimization  process  is  needed  to 
determine  the  appropriate  material  properties  and  thickness  of  the  capsule 
wall  for  maximum  specific  performance  (performance/weight  ratio) .  One 
measure  of  increased  performance  is  the  time  interval  it  takes  from  total 
comminuting  of  Ceramic  to  failure  of  the  capsule  allowing  confined 
comminuted  material  to  extend  projectile  erosion. 

This  work  focuses  on  the  transient  response  of  a  hollow  capsule 
excited  internally  by  a  time  dependent  pressure  pulse.  Two  models  of  the 
capsule  are  developed 

(a)  A  2-D  model  of  4  strips  connected  at  the  edges  forming  a  hollow 

rectangle,  with  finite,  extent  in  the  plane  of  the  cross-section  and 
infinite  in  the  direction  normal  to  it. 

(b)  A  3-D  axisymmetric  model  of  a  short  thin  cylinder  capped  by  thin 

disks  at  the  cylinder's  ends. 

In  the  2-D  model,  the  transfer  matrix  approach  is  utilized  to  relate  state 
vectors  at  the  two  ends  of  a  2-D  strip.  In  the  3-D  axisymmetric  model,  a 
transfer  matrix  is  also  used  for  the  cylinder,  while  each  disk  is  modeled  as 
dynamic  impedances  relating  force  to  displacement  vectors  along  the 
boundary . 

The  internal  pressure  generated  by  the  expanding  comminuted  material 
is  approximated  by  a  hyperbolic  cosine  bell  shaped  distribution  with  origin 
at  the  center  of  the  faceplates  simulating  peak  pressure  at  impact  and 
pressure  decay  remote  from  impact.  This  approximation  is  consistent  with 
both  a  hydro-acoustic  model  of  the  comminuted  material  based  on  the  small 
amplitude  acoustic  equation  and  a  numerical  distribution  determined  by 

"EPIC"  a  general  purpose  "hydrocode" . 

Presented  in  a  forthcoming  paper,  the  acoustic  model  for  internal 
pressure  requires  the  compliance  of  the  capsule  walls.  An  appropriate 
measure  for  compliance  is  dynamic  mobility  of  the  capsule  that  relates 
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internal  pressure  to  normal  acceleration  at  the  boundary.  Dynamic  mobility 
is  derived  from  frequency  response  adopting  modal  method. 

Finally,  a  parametric  analysis  compares  histories  from  the  2-D  and  3-D 
models  with  emphasis  on  peak  stress  and  its  location.  Also,  frequency 
response  of  driving  point  mobility  yields  its  frequency-averaged  radial 
distribution  serving  as  an  input  to  the  acoustic  model  for  internal 
pressure . 

Since  the  3-D  axisymmetric  model  combines  disk  and  cylinder,  a 
literature  survey  focuses  on  each  component  separately  as  analysis  on 
structures  combining  both  is  rare  as  the  preferred  approach  relies  on 
general-purpose  computer  programs. 

A  large  body  of  references  exists  treating  static  and  dynamic  problems 
of  isotropic,  orthotropic,  stiffened  and  rotating  disks.  A  major  portion 
concerns  frequency  response  while  a  smaller  portion  treats  transient 
response.  Soamidas  and  Ganesan  (1991)  analyze  variable  thickness  polar 
orthotropic  disks.  Xiang  et  al .  (1996)  treat  linear  axisymmetric  frequency 

response  with  concentric  stiffeners.  Karunasena  et  al .  (1997)  treat  the 

static  axisymmetric  response  of  a  disk  with  annular  supports  adopting 
Mindlin's  plate  equations.  Liew  and  Yang  (2000)  solve  the  three-dimensional 
free  vibration  problem  of  solid  and  annular  plates  for  symmetric  and 
asymmetric  modes  by  an  approximate  polynomial-Ritz  method.  El-Raheb  and 
Wagner  (1987) ,  (1994)  study  the  axisymmetric  wave  propagation  in  a  disk  from 

impulse  of  short  duration  adopting  a  plate  flexural  model  and  a  3-D 
elasticity  model.  El-Raheb  and  Wagner  (2001)  then  extend  the  axisymmetric 
flexure  model  to  asymmetric  response  from  off-center  impact. 

For  thin  cylinders,  Wang  et  al .  (1997)  utilize  the  Ritz  method  in  the 

modal  response  of  ring-stiffened  cylinders.  Um  and  Kitipornchai  (1998) 
derive  3-D  elasticity  solutions  to  frequency  response  of  open  cylinders.  Lin 
and  Jen  (2003)  adopt  the  Chebyshev  collocation  method  in  the  analysis  of 
laminated  anisotropic  cylinders.  El-Raheb  and  Wagner  (1989)  treat  wave 
propagation  in  a  thin  cylinder  with  concentrated  masses  attached  to  its 
surface  adopting  a  modal  approach.  El-Raheb  and  Wagner  (1985) ,  (1986)  and 

(1989)  adopt  transfer  matrices  to  treat  frequency  response  of  cylinders 
connected  to  toroidal  shells  stiffened  by  discrete  rings  and  annular  disks 
and  damped  by  constrained  visco-elastic  layers  applied  to  the  surface. 

2 .  2-D  model 
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Fig.  2(a)  depicts  the  capsule  geometry  in  the  2-D  model.  In  this 
model,  a  rectangular  plate  with  finite  extent  along  the  local  X  axis  and 
infinite  length  along  the  local  Z  axis  is  termed  a  2-D  strip.  A  schematic 
of  strips  joining  at  corners  through  mobilities  is  shown  in  Fig.  2(b)  where 
reference  coordinate  system  (A,  K)  is  shown  as  dashed  lines.  Corner  l 


(enclosed  in  a  circle)  joins  three  branches  7,7  +  1  and  7  +  2,  through 
linkage  mobilities  ytJ.,  y l]+]  and  y  jJ+2  respectively  where  (i,j)  stands  for  the 

corner-branch  dyad.  In  the  frequency  domain  where  time  dependence  is 
sinusoidal  with  frequency  CO 


z(“}  =(-ma>2  +k+'uog)[j>  ,  a- 1,2,  i  =  V— 1 

and  z,-  ■  is  a  2x2  diagonal  impedance  matrix  with  components  the  one-degree 
of  freedom  oscillators  with  mass  m  ,  spring  k  ,  viscous  damper  g  ,  along  the 


local  coordinates  ( X,y )  of  the  jth  branch.  Let  Ai;  be  the  4x4 
transformation  matrix  relating  local  to  reference  state  vectors  Sloc  —>  S ref 

s5  =  Aus“  >  Su-C-sfc.  Au  = 

where  {f,g  >L  are  force  and  displacement  vectors  at  the  end  of  branch  (i,j) 
and  a,.  •  is  the  2x2  transformation  sub-matrix.  Define  f  and  g  as 

f  =  {fx,qydnz}T  ,  g  =  {ux,uy,y/2f  (3) 


a 


0  a 


0 


(2) 


fx,qv,mz  are  extensional  force,  shear  force  and  moment,  and  ux,uy,y/z  are 
corresponding  displacements  and  rotation  along  local  axes  depicted  in  Fig. 
2(c)  .  Let  /  ,  be  branch  length  then  the  state  vector  at  X  =  lt]  is  related  to 

that  at  X  -  0  by  the  transfer  matrix 


suW- 


+01 

.8(01 


:XJ(o-»OSu(o)- 


111  ^12 

.^21  t22 


(4) 


For  a  2-D  strip  the  transfer  matrix  is  derived  in  Appendix  I.  Since  linkage 
mobilities  y  and  corner  impedances  Z  connect  branch  ends  to  rigid  corners, 


S  at  a  rigid  corner  is  related  to  S  at  a  branch  end  by 
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ff] 

I 

o' 

I 

£e7‘ 

ff]  f 

>  = 

UJ,,; 

i 

Uj 

0 

I 

i>j 

UJ,,;  1 

1  at  x  =  0 
- 1  at  x  =  /  . 


(5) 


At  each  corner,  equilibrium  of  f  and  compatibility  of  g  take  the  form 


iauKj£e=°  ,  i  =  \ne 


;= 1 


(6) 


a i, g u  =  , gM  ,  j  =  2, ,  i'  =  l, nc 


where  is  number  of  branches  joined  at  corner  /  and  Hc  is  total  number  of 

corners.  Substituting  (5)  in  (6)  and  expressing  S in  terms  of  S;  (0) 
using  (4)  produces  6«c  simultaneous  equations  in  the  global  SG(0),  the 
ensemble  of  all  Si ;  (0)  in  the  configuration 

Mg  Sg(0)  =  0  ,  SG(0)  =  {S,  ;(0)}r  ,  (j=\,nbi  ,  i  =  \,nc)  (v) 

Transient  response  is  determined  by  modal  analysis.  The  eigenset  of 
(7)  derives  from 

det[MG]  =  0  =>  {$,0)}*  (8) 

is  the  k‘h  eigenset  formed  of  the  eigenvector 

(8a) 

and  corresponding  eigenvalue  C0k  .  The  modal  analysis  proceeds  by  expanding 
the  global  displacement  vector  gG(x,/)  in  terms  of  the  orthogonal  set  (8a) 

UoM  =  Hbk(t)  ®k(x)  o) 

k 

Substituting  (9)  in  the  forced  elasto-dynamic  equations  of  motion  and 
enforcing  orthogonality  of  the  {«*(*)}  set  yields  uncoupled  ordinary 

differential  equations  in  bk(t ) 

bk  (0  +  co2kbk  (0  =  Na  f(t)/Nkk 

NJk  =  (p(x)\p<bk(x)),  Nkk=(Ok(x)\p®k(x)) 

(')  is  time  derivative,  p(x)  is  spatial  vector  of  external  excitation,  f (t) 
is  its  time  dependence,  and  p  is  density.  Equation  (10)  admits  the  solution 

N  'r 

bk (0  = - ]  f(r)  sin ( 64 (/ -r))dz  (ID 
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3.  3-D  axi symmetric  model 

The  3-D  axisymmetric  model  consists  of  two  thin  disks  capping  the  two 
ends  of  a  short  cylinder  (see  Fig.  3(a)) .  The  cylinder  transfer  matrix  is 
derived  in  Section  3a.  Dynamic  impedance  of  the  disk  at  its  boundary  is 
derived  in  Section  3(b).  Finally,  the  procedure  of  coupling  motions  of 
cylinder  and  disks  is  derived  in  Section  3c. 

3a.  Cylinder  transfer  matrix 

The  cylinder  has  radius  ac  ,  length  lc  ,  plate  thickness  hc  ,  modulus 


E  Poisson  ratio  V  ,  and  mass  density  pc  .  Forces  and  displacements  of  a 

C  1  c 

cylinder  element  are  shown  in  Fig.  3(b) .  For  periodic  motions  in  time  with 
frequency  (0  ,  the  coupled  axisymmetric  elasto-dynamic  equations  of  the 
cylinder  are 


d  N  =  -pho)2ur 

x  xx  r  c  c  c 

dxQx + —  =  -PcKc°2wc 

a c 

Qx  =  dxMxx^~ZO)2dwc  , 


(12! 


X-  1  + 


2k 

1-1/ 


(x,0)  are  axial  and  circumferential  coordinates,  (uc,Wc)  are  axial  and  radial 
displacement,  the  factor  %  multiplying  rotatory  inertia  accounts  for  shear 
deformation  and  K  is  shear  constant  (El-Raheb  and  Wagner  (1998) )  ,  and 
(N  0  M  )  are  axial,  shear  and  moment  resultants  relating  to 
displacements  by  the  constitutive  law 


N, 


EA 


l-V 


o  Wc 

5xuc-vc  — 


N, 


ee 


EA  I  o 

— ^-trl  V  OM, 

l-v2 


a„ 


(13) 


3 

'c§"c 


Eh 


2Z°~WC 


12(1- v2) 

Substituting  (13)  in  (12)  leads  to  a  sixth  order  system  of  ordinary 
differential  equations  in  ( UC,WC ) 


(14) 
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Since  the  operator  Dc  has  constant  coefficients,  the  solution  is  in  terms 
of  exponentials 


W*)J 


CJ 


(15) 


Substituting  (15)  in  (14)  yields  a  homogeneous  set  of  simultaneous  equat: 


Uc 

W 


0  =>  det[Lc]  =  0 

The  implicit  eigenvalue  problem  (16a)  yields  the  eigenvalues  XCj  , 
eigenvectors  [CcuJ,Cc^  ,  and  solutions  for  Uc,Wc 

n,(i)  =  t-Q,eh  - - - - 

;=i  ;=i 


(16a) 


.  (16b) 


Define  the  state  vector  Sc(x)  as 

Sc(x)  =  - 


fc<X>  1 


(17) 


)' 


(18) 


tc  =  {N»>QxxM*}  >  ic  =  K»1 

Substituting  (16b)  in  (13)  then  in  (17) 

Sc(x)  =  Bc(x)  Cc,  Cc  -  {(- cu\X-  cu2’  ■  ’Ccw\>C CW2’  ■■} 

Bc(x)  is  a  matrix  of  exponentials.  Evaluating  (18)  at  the  two  boundaries  of 
the  cylinder  X  =  0 ,lc  then  eliminating  Cc 

Sc(/c)  =  Bc(/c)  B'’(0)-Sc(0)  =  Tc(0  — » lc)  Sc(0)  (19a) 

T  (0  —>lc)  is  the  cylinder  transfer  matrix  relating  Sc  at  its  two 
boundaries.  Eq.  (19a)  can  be  expressed  in  terms  of  four  sub-matrices 


1(4)1 
sll )  I 


tcli  tcU 

^c2l  *c22 


i  1(0)1 

1 8,(0)  I 


(19b) 


3b.  Disk  impedance  matrix 

Each  disk  has  radius  ac  which  is  the  same  as  the  cylinder's,  plate 

thickness  hd  ,  and  modulus  Ed  ,  Poisson  ratio  Vc  and  mass  density  pd  .  Forces 
and  displacements  of  a  disk  element  are  shown  in  Fig.  3(c) .  For  periodic 
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motions  in  time  with  frequency  (O  ,  the  in-plane  axisymmetric  equation  of  the 
disk  is 


dNrr  +  (Nrr  N<^  =  Pdhd®\ 


(20’ 


(r,0)  are  radial  and  circumferential  coordinates,  Ud  is  radial  displacement, 
and  (Nrr ,  N00 )  are  radial  and  circumferential  stress  resultants  relating  to 
displacements  by  the  constitutive  law 

U . 


N  =£A 


N. 


1  ~v] 

EdK 


oe  ,  2 

l-V'd 


dr»d+vd  — 


„  w , 
+  — 


(21) 


Substituting  (21)  ir.  (20)  yields 


[v?  +  pA«2]^  =0> 


(22) 


(23) 


with  solution 

ud(r)  =  CduJl(kder) 

,  _<o_  I  Ed 

Kde~  >  Jn  „2>,  . 

Mindlin' s  axisymmetric  flexural  plate  equations  of  the  disk  are  (El-Raheb 
and  Wagner  (2001) ) 


Dr 


(i-^)vy<f+(i+v<f)a^]-KGA(^+5rw<,)  =  -^t©V</ 


E  h 

«gM ?  V,  +  ) = -MA  • 


(24) 


with  constitutive  relations 


=  D 


drYd  +vd—  >  Qr=  KGdhd  (drwd  +  Yd  ) 


V 


(25) 


Wd,Yd  are  axial  displacement  and  rotation,  Qr,Mrr  are  shear  and  moment 
resultants,  and  K  is  shear  constant.  Eliminating  \)/r  in  (24)  yields  the  4th 
order  equation  in  Wd 
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.2  k 


v2+- 

Y  0  ^  2 

V  C*'J 


. 2  ^ 


V  +  — 
V0  +  2 


\2co2 
c 2  h 2 

Gde'ld 


”7=0,  c*  ,  =. 


* 9l 

Pd 


(26) 


with  general  solution  in  terms  of  Bessel  functions 
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(27) 


^(O  =  cdwj0(V) 


Substituting  (27)  in  (26)  yields  the  dispersion  relation 


(~K  +Ks)(rK  +K)~ 


,  9  12  At; 

2  ,  7.2  \  1  d£_  _  Q 


(28) 


kdc  =  o)/cde,  kds=colcds 

The  four  roots  .  of  (28)  correspond  to  four  primitives  where  only  the  two 
regular  at  r  =  0  are  retained 


wAr)  =  Hc^J  o(V) 

7=1 

¥d 0)  =£Cdv,;  J0  (v) ,  CdV]  =  -y-(^2  +^) 

7=1  ^7 

(  )'  is  derivative  with  respect  to  the  argument. 


Define  the  state  vector  Sd  combining  in  and  out  of  the  plane 
variables 

S„  =  {f,.S,}r  ^  ^  ,30) 

f d={Nn,QrMn}r  >  g d=^„Wd>¥d}T 

Substituting  (23)  and  (29)  in  (21)  and  (25)  yields 

W  =  »*(>■)  C„,  gJW  =  B„(r)C, 

1311 

Cd  =  B gd(cic)  gd(ctc) ,  Cd  =  {CA , Cdwl , Cdw2 } 

and  B  o  are  3x3  matrices.  Eliminating  Cd  from  (31)  determines  the 

Jd  ga 

disk  impedance  matrix 

=  g(r)  =  ZJri  g  /r)  132) 

is  the  impedance  matrix  of  the  disk  at  the  boundary . 


3c..  Coupling  cylinder  and  disks 

Continuity  of  S  at  the  cylinder-disk  interface  X-Xe=0,lc  and  r  =  Clc 
requires  expressing  gd}  in  the  coordinate  system  of  {fc>gc}  bY  means  of  a 
re-ordering  matrix 


12  ■ 


(33) 


fe(x.)  =  Rff d(ac), 


Rf 


0-10 

1  0  0 

0  0-1 


gc(^)  =  -Rf  %d(ac)  =>  gd(a'C)  =  -Rf1gC(XJ 


Substituting  (32)  in  (33)  then  eliminating 

t(o)=z„fc(o),  f,(y  =  -z,g,(t)  |34| 

Z,  =  -Rf  Z„K)  Rf1 

Zf .  is  disk  impedance  in  the  cylinder  coordinate  system.  Applying  (34)  at 


the  two  boundaries  of  the  cylinder  and  relating  fc(4),gc(4)  to  fc(0) >  gc(0) 
using  (19b)  produces  the  homogeneous  equations 


1 

' 

0 

0 

M  = 

tel. 

tc12 

-I 

0 

(35) 

c 

^c2\ 

^c22 

0 

-I 

o 

o 

I 

zj 

and  O  are 

3x3 

unit 

and  null  matrices, 

and  0 

Mc  is  a  12x12  matrix, 
the  null  vector  of  order  6.  A  non-trivial  solution  of  (35)  yields  the 
implicit  eigenproblem 


det[Mc]  =  0  =>  {$„»}, 


(36a) 


Eq.  (31)  and  (33)  determine  the  disk  eigenfunction  <E>d(r)  .  The  global 
eigenfunction  of  cylinder  and  conjoined  disks  then  follows 

=  l36b> 

S  is  an  intrinsic  coordinate  along  the  boundary;  S  =  X  on  the  cylinder  and 
j's/'  on  the  disks.  Modal  analysis  proceeds  by  expanding  the  global 

displacement  vector  gG(^,0  =  {g,(4  g in  terms  of  the  orthogonal  set 


(36b) 

&,(*,') =EA*w®«w  1371 

k 

Substituting  (37)  in  the  forced  elasto-dynamic  equations  of  motion  of 
cylinder  and  disks  and  enforcing  orthogonality  of  the  { ^ok (s) /  set  yields 

uncoupled  ordinary  differential  equations  in  bk(t) 
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(38) 


bk(t)  +  co2kbk(t)  =  Nfkf(t)/Nkk 
Nfk  =  <pCv)  |  p%k(s)) ,  Nkk  =  <®ot(j)  |  pOGi0)> 

(')  is  time  derivative,  p(s)  is  spatial  vector  of  external  excitation,  f (t) 
is  its  time  dependence,  and  p  is  density.  Equation  (38)  admits  the  solution 


K(  0  =  - 


N 


fk 


k  Nkk  0 


f  /0)  dz 


(39) 


For  frequency  response,  bk(f)  —  —  CO1  bkQ  reducing  (38)  to 

K(<)  =  N,le“l(Nkl(-<o1+al))  (40) 

CO  is  frequency  of  excitation. 


4 .  Forcing  pressure 

In  the  3-D  axisymmetric  model,  a  spatial  distribution  of  the  applied 
pressure  is  approximated  by  the  bell-shaped  function 

p(r)=  p0  +  P\  sech(p2r)  ,  0  <r<ac  (41) 

where  the  triad  (p0 ,p},p2)  of  constants  characterizes  translation, 

magnification  and  maximum  slope  of  the  profile.  Different  triads  are 
assigned  to  top  and  bottom  disks  (see  Fig.  3(a)) .  A  typical  normalized 
distribution  is  shown  in  Fig.  4. 

In  the  2-D  model,  assuming  that  pressure  acts  over  a  distance  lx  along 
Z  where  /  is  length  of  capsule  along  top  or  bottom  faces  (see  Fig.  2(a) ) ,  a 
profile  similar  to  (41)  is  used  but  with  Pq,  P\  scaled  to  conserve  total 
force 

P(x)=  Po  +  p\  sech(p2(x  -0.5  lx)) ,  0  <  x  <  lx  (42) 

The  condition  of  force  conservation  yields 
lx  ac 

lx  J"  p(x)  dx  -  2n  J  p{r )  rdr  ,  lx=  2 ac  (43) 

0  0 

Assuming  that  p()  /  p^  =  p0  /  /">,  ,  then  Eq .  (43)  suffices  to  determine  the  dyad 

(Po>Pi)  ■ 
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6 .  Results 


The  2-D  capsule  has  a  rectangular  cross-section  made  of  4  plates  as 
shown  in  Fig .  5 : 

(1)  top  plate  joining  corners  1  through  5 

(2)  bottom  plate  joining  corners  6  through  10 

(3)  left  plate  joining  corners  1  and  10 

(4)  right  plate  joining  corners  5  and  6 

A  global  coordinate  system  is  defined  for  both  the  2-D  and  3-D  models  with 

origin  at  corner  1  such  that  {X:  0  <  X  <  /^j  0  0  <  Y  <  Wy  j  where  (lx ,  Wy )  are 

length  and  width  of  the  capsule.  The  following  geometry  and  material 
properties  are  considered 

lx  =  100  mm  ,  w  =  12.5  mm  ,  h  -  2.54  mm 
E  =  2  xlO12  dyn/cm2 ,  p  =  %g/cm3,  h  =  0.3 

where  h  is  uniform  plate  thickness  and  ( E,p,V )  are  Young's  modulus,  density 
and  Poisson  ratio.  In  the  2-D  model,  3  stresses  are  monitored  along  each 
segment's  local  X-axis:  flexural  stress  O xxJ-  ,  shear  stress  Tx2  and 

extensional  axial  stress  O xxe  .  In  the  3-D  model  4  stresses  are  monitored; 

either  0’xxf,Tja,<Jxxe,(Tsge  along  the  cylinders  X-axis  which  is  the  global  Y- 

axis,  where  Ggee  is  the  extensional  circumferential  stress,  or  else 

<Jrrf ,  Xn ,  <7ne ,  Gg&e  along  the  disk's  radius  which  is  the  global  X-axis.  Stress 

at  a  point  with  global  coordinates  (X,F)  will  be  referred  to  as  (j{X,Y)  . 

The  first  10  mode  shapes  and  corresponding  resonant  frequencies  are 
shown  in  Figures  6(a-j).  The  fundamental  at  1.063  KHz  (Fig.  6(a))  is 
symmetric  about  the  axis  of  bilateral  symmetry.  This  corresponds  to  the 
first  symmetric  flexural  mode  of  the  two  faces.  The  second  mode  (Fig.  6(b)) 
is  the  first  anti - symmetric  flexural  mode  of  the  faces.  Higher  flexural 
modes  with  larger  wave  numbers  are  shown  in  Fig.  6(c-f)  and  then  the 
fundamental  symmetric  extensional  mode  at  21.72  KHz  (Fig.  6(g)),  followed  by 
the  first  anti -symmetric  extensional  mode  (Fig.  6(h))  at  24.54  KHz.  At 
higher  frequencies,  (Fig.  6  ( i - j ) )  the  flexural  modes  exhibit  noticeable 
motion  of  the  side  plates,  which  were  evanescent  for  frequencies  below  the 
fundamental  symmetric  extensional  mode  in  Fig.  6(g) .  For  each  plate,  modes 
belong  to  one  of  4  groups  depending  on  symmetry  of  the  mode  (symmetric  or 
anti -  symmetric)  and  type  of  mode  (flexural  or  extensional)  . 
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In  the  3-D  axisymmetric  model,  a  spatial  distribution  of  the  applied 
pressure  on  top  and  bottom  plates  is  approximated  by  the  bell-shaped 

function  in  Eq .  (38)  where  the  triad  of  constants  ^Pq->P\^PP)  characterizes 

translation,  magnification  and  maximum  normalized  slope  of  the  profile. 
Different  triads  are  assigned  to  top  and  bottom  plates  (see  Fig.  4) .  The 
normalized  peak  pressure  for  the  axisymmetric  excitation  is  unity.  In  the  2- 

D  model,  the  equivalent  bell  shaped  distribution  with  (p0,px)  that  conserves 

total  force  is  determined  applying  Eq .  (39)  .  The  resulting  ( Po>P\)  i-s  listed 

in  Table  I  showing  that  in  the  2-D  model,  peak  applied  pressure  on  the  top 

plate  is  p  =  p0  +  P\  ~  0  4  •  The  top  plate  is  divided  into  4  equidistant 

intervals  of  12 . 5mm  along  X  while  the  side  plate  is  divided  into  2 

equidistant  intervals  of  6.25mm  along  Y  .  The  time  dependence  of  the 
pressure  pulse  is  trapezoidal  lasting  25  /US  with  5  jUS  rise  and  fall  times 

and  a  15  JUS  plateau. 

Fig.  7  plots  stress  histories  on  the  left  plate  for  0  <Y  <12.5  mm  .  The 
peak  of  a ^ (0, 0)  reaches  38  pm  (Fig.  7(a))  pointing  to  a  location  where 

material  may  yield  and  ultimately  fail.  <J  ^  (Fig.  7(c))  is  relatively  small 

since  it  reacts  to  the  difference  between  integrated  pressure  over  the  top 

plate  («  /x /4  =  (9(1) )  and  its  inertia.  Fig.  8  plots  stress  histories  on  the 

top  plate  for  12.5  <X  <50  mm  at  12.5mm  intervals.  The  peak  of  0^(50, 0) 

(Fig.  8(a))  reaches  23  pmx  pointing  to  another  location  of  considerable 

flexural  stress.  Fig.  8(b)  plots  histories  of  shear  stress  X x2  with  peaks 

not  exceeding  pnx  .  Fig.  8(c)  plots  histories  of  extensional  stress  <7 ^  in 

the  plane  of  the  plate  with  a  peak  of  CT^ (50,0)  reaching  5  .  It  is  evident 

that  the  two  most  vulnerable  locations  with  high  <7^  are  at  the  interface 

of  top  and  side  plates  and  at  the  top  plate's  center.  Extensional  and  shear 
stresses  do  not  contribute  to  yielding  and  failure  of  the  material  unless 
deformation  is  sufficient  to  promote  additional  stress  from  nonlinear 
stretching . 

Geometry  and  material  properties  in  the  3-D  model  replicate  those  in 
the  2-D  model.  The  first  6  mode  shapes  are  plotted  in  the  left  column  of 
Fig.  9  while  the  corresponding  modes  from  the  2-D  model  are  plotted  in  the 


16 


right  column  for  comparison.  The  fundamental  mode  1  at  2.075  KHz  is  the 
first  symmetric  flexural  mode  of  the  disks  while  mode  2  at  4.714  KHz  is  the 
first  anti - symmetric .  Flexural  modes  3,4  and  6  then  occur  before  the 
fundamental  symmetric  extensional  mode  5  occurs  at  29.74  KHz.  The  frequency 

ratio  C03_D  /  CO^_D  for  modes  1,2  and  4  is  approximately  2  (Fig.  9)  .  This  ratio 

drops  to  less  than  1.3  for  all  following  modes  suggesting  that  the 
difference  in  response  between  3-D  and  2-D  models  is  largest  at  low 
frequencies.  It  also  implies  that  for  a  pulse  of  short  duration  compared  to 
the  fundamental  period,  transient  response  from  the  two  models  may  be 
comparable  since  then  response  would  be  controlled  by  modes  with  short  wave 
length . 

Fig.  10  plots  stress  histories  on  the  cylinder  for  0  <Y  <12.5/72/72  .  The 
peak  of  0^(0, 6)  (Fig.  10(a))  is  half  that  in  2-D  while  the  peak  of  O ^ 

(Fig.  10(c))  is  almost  the  same  as  that  in  2-D.  0’eSe  (fig.  10(d))  is 
approximately  l/4Cr^  and  has  no  counterpart  in  2-D.  Fig.  11  plots  stress 

histories  on  the  top  disk  for  0  <  X  <  50/77/7/  at  12.5mm  intervals.  The  peak  of 
O'  rrj  (50, 0)  reaches  the  same  magnitude  as  that  of  0*^(0, 0)  in  2-D,  and  the 

peak  of  Orrf(0, 0)  reaches  the  same  magnitude  as  that  of  <V(50>°)  in  2-D. 

Fig.  11(b)  plots  histories  of  shear  stress  T x2  with  a  peak  1/3  that  in  2-D. 
a rre  and  Ogge  (Fig.  11(c),  (d)  )  have  approximately  the  same  magnitude  reaching 

1/5  that  of  O ae  in  2-D.  Consequently  in  3-D,  magnitude  of  peak  flexural 

stress  is  approximately  the  same  as  in  2-D  but  locations  are  reversed; 
highest  stress  is  at  the  disk's  center  followed  by  that  at  the  edge.  It  is 
surprising  that  in  spite  of  the  different  dimensionality,  histories  of 
flexural  stress  from  the  two  models  agree  qualitatively  and  quantitatively. 
The  discrepancy  between  the  two  models  is  in  the  extensional  stress  that  is 
substantially  lower  than  flexural  stress. 

Fig.  12  follows  the  evolution  of  deformation  of  the  3-D  capsule  for 
10  <t  <80/75.  The  dominant  mode  of  deformation  is  a  symmetric  bulging  that  is 
larger  on  the  top  plate  than  on  the  bottom  due  to  the  difference  in 
magnitude  of  the  applied  pressure. 
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7 .  Conclusion 

Two  models  are  presented  for  transient  response  of  a  thin  walled 
capsule  forced  by  a  time  dependent  internal  pressure.  The  2-D  model  has 
finite  dimensions  in  the  plane  of  the  cross-section  and  extends  to  infinity 
normal  to  this  plane.  The  3-D  axisymmetric  model  includes  a  short  thin 
cylinder  capped  by  two  thin  disks.  Both  models  adopt  transfer  matrices 
relating  state  variables  at  two  ends  of  a  segment  and  transient  response 
based  on  modal  analysis.  Noteworthy  results  are 

1)  For  low  wave  numbers,  eigen- frequencies  of  the  3-D  model  are  twice 
those  of  the  2-D  model.  The  difference  diminishes  as  wave  number 
rises  . 

2)  The  assumed  bell-shaped  internal  pressure  distribution  agrees  with 
that  determined  from  an  acoustic  model  and  from  a  numerical  general 
purpose  program. 

3)  In  both  the  2-D  and  3-D  models,  peak  flexural  stress  occurs  at  the 
junction  of  cylinder  and  disks  as  well  as  center  of  disks. 

4)  Peak  extensional  and  shear  stresses  are  one  order  of  magnitude  smalle 
than  peak  flexural  stress  except  for  circumferential  stress  in  the 
cylinder  with  peak  1/4  that  of  flexural  stress. 

5)  In  spite  of  the  different  dimensionality,  histories  from  the  two 
models  are  in  good  agreement. 


Appendix  A.  Transfer  matrix  of  2-D  strip 

For  sinusoidal  time  dependence  with  frequency  CO  ,  the  linear  elasto- 
dynamic  equations  of  a  2-D  strip  are 


d2us 

dx2 


~k2ssus  ,  kse  =  <o! cse  ,  cse=^Es/ps( \-v2) 


D 


dx' 

KG'h 


-kGA 


Ys+- 


V 


dw_ , 

dx 


^——co:y/ 

12 


(Al) 


f82ws  |  djQ 

dx2  dx 


=  ~aA®2w  ,  1) 


EM 


u(i -k) 


11  ^ ,  \Vs ,  y/s  are  displacements  along  local  coordinates  X  and  y  (see  Fig.  lc 


E,,GS,V,,PS  are  extensional  and  shear  moduli,  Poisson  ratio,  and  mass 

density,  h 5  is  thickness  and  K  is  shear  constant.  The  corresponding 
constitutive  relations  are 
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/* 


Eshs  dus 

O-^2) 


dw 


cL  =  kGA\  -t f  +  ¥s  ,  mx  =  D 

ox  y 


dWs 

dx 


( A2 ) 


/  , qx , mx  are  axial  force,  shear  force  and  moment  (see  Fig.  lc)  .  The  general 


solution  of 

(Al) 

< 

w,(*) 

»  = 

k,rx 


0 


0 

A.yX 


0  S,tel-  S,2e 


Ayl  2  satisfy  the  dispersion  relation 

>2 


\c  1 

< 

Gsw\ 

c  , 

l  sw2  J 

,  S>j=-y(K+kl)  <A3> 


_ _ _ 1? ir2  _ 

(-F  +Yr)(-Y.  +  k2 J-  — ~  =  6 

\  SJ  S£S\  SJ  sss  ^2 


(A4) 


( A5 ) 


K  =  >  c„=4kGJp* 

Substituting  (A4)  in  (A2)  in  terms  of  the  state  vector  Sf  =  {fs. ,  g,v } 

S,(x)  =  B,(x)C, 

f,  ={fx,qx>™x}Tt  , 

B,.(x)  is  a  3x3  matrix  of  exponentials  and  Cs  —  {C4.a  ,  CwJ ,  Csw2  }  ■  Evaluating 

(A5)  at  the  two  ends  of  the  strip  (0,4)  then  eliminating  Cs  yields  the 
strip' s  transfer  matrix 


S,  (4 )  =  Bs  (4 )  b;1  (0)  s,  (0)  =  ts  (<>-►/,)  s.  (0) 


( A6 ) 
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Table  I.  Parameters  of  pressure  profile 
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4.  Radial  distrib ; 


top  plate 


Figure  5.  Cross-section  of  2-D  capsule  with  comer  and  branch  numbers 


(c)  0=  5.941  kHz  (h)  n=  24  545  kHz 


(d)  n=  7.627  kHz  0  27-252  kfiz 


(e)  n=  14.737  kHz  (j)  0=  30.596  kHz 


Figure  6.  Resonant  frequencies  and  modes  of2-D  capsule 


3-D  axisymmetric 


2-D 


n=  2.075  kHz  made  1  0=  1.063  kHz 


0=  4.174  kHz  mode  2  0=  1.907  kHz 


fl=  8. 317  kHz  mode  3  Q=  5.941  kHz 


n=  14.572  kHz  mode  4  0=  7.627  kHz 


0=  18.762  kHz  mode  5  71=  14.737  kHz 


n=  29.741  kHz  mode  6  Q  =  21.724  kHz 


Figure  9.  Comparison  of  modes  from  3-D  axisymmetric  and  2-D  models 


t  (\is)  t  (\Ll) 

Figure  10.  3-D  stress  histories  on  cylinder:  (a)  axxf,  (b)  xx2 ,  (c)  axxe ,  (d)  a00e 
- Y=0,  ------ Y  =  0.25  wy,  .  Y=0.5wy, 

- Y  =  0.75wv,  -  -  -  Y  =  wv 

y  y 


Figure  11.  3-D  stress  histories  on  top  disk:  (a)  arrf,  (b)  xxr ,  (c)  arn 


- X-  0.5  lx,  . X  =  0.375  lx,  - X  -0.25  l} 

-  -  -  X  =0.125  L,  -  -  -  X=0 


